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Kyoto University, Kyoto 606-8502, Japan 

We review recent developments in the theory of inflation and cosmological perturbations 
produced from inflation. After a brief introduction of the standard, single-field slow-roll 
inflation, and the curvature and tensor perturbations produced from it, we discuss possible 
sources of nonlinear, non-Gaussian perturbations in other models of inflation. Then we de- 
scribe the so-called 8N formalism, which is a powerful tool for evaluating nonlinear curvature 
perturbations on super Hubble scales. 



I. INTRODUCTION 

One of the most successful applications of the theory of general relativity is cosmology. Over the 
past half century the big-bang theory of the universe, that the universe was born in an extremely 
hot and dense state, expanded explosively and cooled down to the present state, was observationally 
tested from various aspects and it is now firmly established. According to the big-bang theory, our 
universe is about 14 Giga years old, and the universe was radiation-dominated in the beginning. 
It became matter-dominated when the universe was about 100,000 years old, which happens to be 
about the same time when the photons decoupled from baryons, and started to travel freely until 
today, which are observed as the cosmic microwave background (CMB) radiation. The epoch when 
the CMB photons were scattered last before they reach us forms a 3-dimensional hypersurface, and 
it is called the last scattering surface (LSS). 

Despite its tremendous success, there are still a couple of very basic problems that the big-bang 
theory cannot explain. One of them is the horizon problem or perhaps better to be called the 
causality problem, and the other the flatness problem or the entropy problem. 
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A. Horizon problem 

Let us first consider the horizon problem. The big-bang theory assumes an homogeneous and 
isotropic universe on large scales. So the metric is assumed to be in the form, 

ds 2 = -dt 2 + a 2 (t)daf 3) , (1) 



where da 2 3 s is the 3-metric of a constant curvature space with K being the curvature, ^R^km = 
K(S t k 5m — 5^5^). A coordinate system that spans dcr^) ^ s sa id to be comoving because an observer 
staying at a fixed point on the 3-space is comoving with the expansion of the universe. In this 
spacetime, the time-time component of the Einstein equations, the Friedmann equation, is 

where M 2 X = (87rG) _1 in the units H = c = 1, and the trace of the space-space components of the 
Einstein equations gives 

a p + 3P 



a~ 3M 2 pl ' 



(3) 



where p is the energy density and P is the pressure of the universe. This latter equation shows 
that the expansion of the universe is always decelerating as long as p + 3P > 0, which holds for 
both radiation P = p/3 and matter P = 0. For simplicity, if we assume a simple equation of state 
P/p = w =constant and K = (which should be a good approximation in the early universe when 
w = 1/3 since p oc a~ 3 ( 1+w ^ = a~ 4 ), one finds 

aocf; n =w ^<l for w > -I . (4) 

This result may be regarded as a consequence of the attractive nature of the gravitational force. 



Now we introduce the conformal time dr\ = dt/a(t), and rewrite the metric as 

ds 2 = a 2 (r])ds 2 ; ds 2 = -drj 2 + da 2 3) . (5) 

Since the conformal transformation of the metric does not change the causal structure, the static 
metric ds 2 perfectly describes the causal structure of the universe. If the range of rj were infinite 
to the past, there would be no horizon problem. The problem is that the conformal time is finite 
in the past if w > — 1/3 or p + 3P > 0, because 

r dt ' f tdt ' 2 ^ 
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FIG. 1: Horizon problem. The conformal time of the last scattering surface t]lss from 77 = is about 1/30 
of that of today 770- 



This implies that the size of lightcone emanating from a point at the beginning of the universe 
when f] = will cover only a finite fraction of spacetime. Since the comoving distance traveled 
by light is equal to the corresponding conformal time interval, the comoving radius of the causally 
connected region on the LSS is equal to its conformal time t]lss- From the fact that the LSS is 
located at redshift z ~ 10 3 and the universe is approximately matter-dominated since then, one 
finds that this region will cover only a tiny fraction (about 10~ 3 sr) of the sky. This is the horizon 
problem (see Fig. 1). 

The solution is clear: The horizon problem disappears if the conformal time is either infinite in 
the past or the beginning of the universe r\ = is extended sufficiently back in time to cover the 
whole visible universe. Since the comoving radius of the visible universe on the LSS is r/o — r/LSS 1 
where r/o is the conformal time today, the problem is solved if tjlss > Vo ~ Vlss- I n Einstein 
gravity, this means that the equation of state must be w < —1/3 or the expansion of the universe 
must be accelerating (a > 0) for a sufficient lapse of time in the very early universe. 

Here we should note that solving the horizon problem does not mean explaining the homogeneity 
and isotropy of the universe. As clear from the above argument, we had to assume the homogeneity 
and isotropy of the universe to pose the horizon problem. This point is very often misunderstood 
in the literature. 



B. Flatness problem 

Again we assume a spatially homogeneous and isotropic universe, Eq. (1). The Friedmann 
equation (2) tells us that the curvature term K/a? is completely negligible in the early universe 
when p oc a~ 4 . Conversely, if the curvature term was of the same order of magnitude as the density 
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at an epoch in the early universe, the universe must have either collapsed (if K > 0) or become 
completely empty (if K < 0) by now. 

Alternatively, since the energy density is dominated by radiation in the early universe and so is 
the entropy of the universe, the problem may be rephrased as the existence of huge entropy within 
the curvature radius of the universe, 

S = T 3 (-±=^ 3 « T 3 (-^=) 3 > T$Hl ps 10 87 , (7) 

where To ps 2.7 K is the CMB temperature today [1] and Hq ps 72km/s/Mpc is the Hubble 
constant [2]. Hence the flatness problem may be called the entropy problem. 

It is then apparent that the solution to the flatness problem needs huge entropy production at 
a sufficiently early stage of the universe. 



C. Inflation as a solution to horizon and flatness problems 

A simple and perhaps the best solution to the horizon and flatness problems is given by the infla- 
tionary universe [3, 4]. Let us assume that the universe was dominated by a spatially homogeneous 
scalar field. For a minimally coupled canonical scalar field (f>, we have 

p=\^ + V{ff>), P=^ 2 -V(cf ) ), (8) 

so p+3P = 2 {(p 2 — V (</>)). Hence ii <ft 2 < V((p), we may have accelerated expansion. In particular, if 
the energy density is dominated by the potential energy, cp 2 <C V ((/>), the motion of the scalar field 
can be ignored within a few expansion times ~ H -1 , and the universe expands almost exponentially, 

H 2 ps — ^-rr ps constant. (9) 
3M 2 pl 

The curvature term K/a 2 becomes completely negligible. 

Thus if the universe is dominated by the potential energy, or the vacuum energy, and the 
potential energy is converted to radiation after a sufficient lapse of time of such a stage, a huge 
entropy is produced and the horizon and flatness problems are solved simultaneously. 



II. SLOW-ROLL INFLATION AND VACUUM FLUCTUATIONS 

There have been a number of proposals for inflationary models. Among others, a simplest class 
of models, and which explains the observational data almost perfectly, is the slow-roll infation [5-7]. 
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The field equation for (f> and the Friedmann equation are 



+ 3H<j) + V'((j)) = 0, H 2 



-6 2 + V(cP) 



where we have justifiably neglected the curvature term. 



(10) 



logL 




<t> = <h 
inflation 
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FIG. 2: The Hubble radius!/ = H 1 and the length scale L — a/k of a comoving wavenumber k in the 
inflationary cosmology, and the definition of the number of e-folds N(<f>). 



The standard slow-roll condition consists of two assumptions. One is that 4> is negligible com- 
pared to 3H(p in the field equation, that is, the equation of motion is friction-dominated. The 
other is that the kinetic term <j) 2 /2 is negligible compared to the potential term V in the energy 
density. Under this condition we have 



i— m 

Y 3H 

Then the potential energy dominance implies 

H ^ 



H 



V 



(11) 



H 2 li 2 



2<r + v 



2V ~ 2 V 2 



ev < 1, 



(12) 



that is, the universe is expanding almost exponentially, and the friction-dominated equation of 
motion \cf>/(3H<p)\ <C 1 implies 



V" 



V" 



3H 2 ^-= W «1 



(13) 



The single-field slow-roll inflation satisfies these conditions. 
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The important property of slow-roll inflation is that Eq. (11) is completely integrable since H 
is a function of 4>. In particular, there is one-to-one correspondence between <j> and t. So instead 
of the cosmic time t we may measure the time in terms of the value of the scalar field. 

Here we introduce a quantity which plays a very important role in the dynamics of slow-roll 
inflation, namely the number of e-folds counted backward in time, say from the end of inflation to 
an epoch during inflation, 



a(t, 



end, 



= exp[N(t -> t end )] -> N = N(<j>) = Hdt. (14) 

Jt(d>) 



a W Jt(4>) 
Its important property is that by definition it does not depend on how and when the inflation 
began. As shown in Fig. 2, N is uniquely determined in terms of the value of the scalar field (up 
to a constant which depends on the choice of an epoch from which N is computed), and one can 
associate iV with the time at which a given comoving wavenumber k crossed the Hubble radius, 
k = aH, at which the value of the scalar field was <fe N = N((/>k). As we shall see below, this turns 
out to be an essential quantity for the evaluation of the curvature perturbation from inflation. 

A. Curvature perturbation 

Let us now consider the curvature perturbation produced from inflation. It arises from the 
quantum vacuum fluctuations of the inflaton field 4>. Since a rigorous derivation would take too 
much space, here we give an intuitive, rather hand-waving derivation. We caution that it could 
well lead to an incorrect result if used blindly. 

The vacuum fluctuations of the inflaton field with a comoving wave number k is given simply 
by its positive frequency function, ipk- Because of the condition V"/H 2 <C 1, on scales k/a 3> H, 
the inflaton field fluctuation behaves like a minimally coupled massless scalar. Hence we have 

(6<i>\k) 2 = \<p k \ 2 , tp k ~ . 1 e-* 1 **; w fe = -»iT. (15) 

As the universe expands the physical wavenumber decreases exponentially and becomes smaller 
than the Hubble parameter, k/a < H, or the physical wavelength exceed the Hubble radius. Then 
the oscillations of ipu are frozen. This could be regarded as "classicalization' of the quantum 
fluctuations. Not that this is merely an interpretation. In a more rigorous sense, freezing of the 
mode function is a process toward infinite squeezing of the vacuum state. 
Setting a = k/H in Eq. (15) gives 
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Therefore the mean square amplitude in unit logarithmic interval of k is 



k/a=H 



(17) 



Inclusion of the non-trivial evolution of the background spacetime and the coupling of the scalar 
field fluctuation with the metric fluctuation do not change the above estimate if we interpret 8(f) 
in the above as those evaluated on the flat slicing, that is, on hypersurfaces on which the spatial 
scalar curvature remains unperturbed. 

It is known that the curvature perturbation on the comoving hypersurface 1Z C is conserved if 
the perturbation is adiabatic [8]. The comoving hypersurface is defined as a surface of uniform 
4>. Then the gauge transformation from the flat slicing to the comoving slicing gives the relation 
between 1Z C and 6(j>, 



n c = — ^5cj) . 



(18) 



Since this is conserved for k/a < H, the spectrum of the comoving curvature perturbation in unit 
logarithmic interval of k is given by 



) k/a=H 



(19) 



A rigorous, first-principle derivation of the above result was first done in [9, 10]. 

The important relation of the above result with the number of e-folds was first pointed out 
in [11]: If we rewrite Eq. (14) as 



N 



^end 



Hdt 



end J-£ 



we find 



5N(<f> k ) 



dN 



H 



(20) 



(21) 



k/a=H 



k/a=H 

provided that we identify 6(f) with the scalar field fluctuation evaluated on the flat hypersurface. 
This is called the 6N formula. 

The 6N formula implies that we only need the knowledge of the background evolution to obtain 
the power spectrum of the comoving curvature perturbation, once we know the amplitude of the 
quantum fluctuations of the scalar field at the horizon crossing (i.e. when k/a = H). It is quite 
generally given by H/ (2ir) in slow-roll inflation. With careful geometrical considerations, the 6N 
formula can be extended to general multi-field inflation [12], 



Vn(k) 



dN ON 



2.) " VAr||2; l' ViV " 2sG °^W 



(22) 
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where G ab is the field space metric and it is assumed that the vacuum expectation values are given 
by 

{5m b )=G ab (f\ . (23) 



2ir r 

The nonlinear generalization of the SN formalism will be discussed in Sec. IV. 

B. Tensor perturbation 

There are not only vacuum fluctuations of the inflaton field but also those of the transverse- 
traceless part of the metric, d l hj^ = S^hf? = 0, that is, the tensor perturbation or gravitational 
wave degrees of freedom. If we construct the second-order action for hfj, we find 

M 2 r 

S-^j d*x^-g{hJT? + -... (24) 

To quantize hj? it is convenient to normalize the kinetic term to the canonical form. This gives 

S ~ \ J d^xV^dij) 2 + • • • ; ^ es ^ T (25) 

If one writes down the field equation for cpij, one finds its mode function (p^ obeys exactly the same 
equation as the one for a minimally coupled massless scalar field, 

k 2 

fc + 3^ fc + — fc = O. (26) 
cr 

Since there are two independent degrees of freedom in fy, the power spectrum of the tensor 
perturbation hj? is obtained as 

m „s 4 4 ^ i , ,2 8F 2 

P ^> = Mr 2 *(2^' 2 =(2W <27) 
Taking the ratio of the tensor spectrum to the curvature perturbation spectrum, we find [12] 

r 4- 8W = "V' (28) 
where tit is the tensor spectral index, ut = dhiVxik) / 'din A;, and the equality holds for the case 
of single-field slow-roll inflation. This is a consistency relation in general slow-roll inflation. As a 
proto-type example, if we consider chaotic inflation [7], we expect to have r ~ 0.1. 

The important point to be kept in mind is that the existence of the vacuum fluctuations of the 
tensor part of the metic is a proof of the existence of quantum gravity. These fluctuations exist in 
any theory of gravity that respects general covariance, apart from possible inessential modifications 
of the spectrum. Thus a clear detection of the tensor spectrum will be a confirmation of not only 
the inflationary universe but also quantum gravity. 
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III. ORIGIN OF NON-GAUSSIANITY 

The standard, single-field slow-roll inflation predicts that the curvature perturbation is a Gaus- 
sian random field and it has an almost scale-invariant spectrum. This seems to fit the current 
observational data quite well [13], it is quite possible that the actual model turns out to be non- 
standard. Maybe it is multi-field, maybe non-slow-roll and/or non-canonical. In such a case, the 
curvature perturbation may become non-Gaussian. Search for possible non-Gaussian signatures in 
the primordial curvature perturbation has become one of the important directions in observation 
in recent years [14]. 

Here we consider possible origins of non-Gaussianity in the curvature perturbation. Essentially 
one can classify the origins into three categories: (1) Self-interactions of the inflaton field and/or 
non-trivial vacua, (2) multi-field dynamics, and (3) nonlinearity in gravity. 

The non-Gaussianities of the first category are generated on subhorizon scales during inflation, 
hence they are of quantum field theoretical origin. Those of the second category are usually 
generated on superhorizon scales either during or after inflation, and they are due to nonlinear 
coupling of the scalar field to gravity. Since they are generated on superhorizon scales, they are 
of classical origin. Finally those of the third category are due to nonlinear dynamics in general 
relativity. Hence they are generated after the scale of interest re-enters the Hubble horizon. Since 
the last category is not really primordial in nature, let us focus on the first two categories. 

A. Non-Gaussianity from self-interaction/non-trivial vacuum 

It is known that conventional self-interactions by the potential are ineffective [15]. This can be 
seen by considering chaotic inflation, for example. In the simplest case of a quadratic potential, 
V = m 2 4> 2 /2, the inflaton is actually a free field apart from the interaction through gravitation 
perturbations. But the gravitational interaction is Planck-suppressed, i.e., it is always suppressed 
by a factor 0{M~ 2 ). In the case of a quartic potential, V = A0 4 , it is known that A should be 
extremely small A ~ 10~ 15 in order for it to be consistent with observation. 

Thus some kind of unconventional self-interaction is necessary. A popular example is the case of 
a scalar field with a non-canonical kinetic term such as DBI inflation [16]. In this case the kinetic 
term takes the form, 

k ~ rHtfyJi-mft = r V 1 • (29) 
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If we expand this perturbatively, 



K = K + S X K + 5 2 K + 6 3 K + ■ ■ ■ , 



(30) 



we will find 



5 2 K oc 7 3 , S 3 K oc 7 : 



,3+2 



(31) 



since 5j = j 3 5X where X = f<j) 2 /2. If we regard the third order part as the interaction, the above 
implies that the scalar field fluctuation will be expressed qualitatively as 



5<f) ~ 5<j)o + 7 2 50o H 



(32) 



where 54>q is the free, Gaussian fluctuation. Thus the non-Gaussianity in 5<p may become large if 
7, which mimics the Lorentz factor, is large [17]. 

A non-trivial vacuum state is another source of non-Gaussianity. If the universe were a pure 
de Sitter spacetime, gravitational interaction would be totally negligible in vacuum, except for 
the effect due to graviton (tensor mode) loops. This may be regarded as due to the maximally 
symmetric nature of the de Sitter space, 50(4,1), which has the same number of degrees of 
symmetry as the Poincare (Minkowski) symmetry. In slow-roll inflation, the de Sitter symmetry is 
slightly broken. Nevertheless the effect induced by this symmetry breaking is small because it is 
suppressed by the slow-roll parameter e = —H/H 2 . 

However, if the vacuum state does not respect the de Sitter symmetry, there can be a large non- 
Gaussianity. Such a deviation from the quasi-de Sitter vacuum, usually called the Bunch-Davies 
vacuum, may occur in various situations, studied e.g. in [18, 19]. 



Non-Gaussianity may appear if the energy momentum tensor depends nonlinearly on the scalar 
field even if the fluctuation of the scalar filed itself is Gaussian. This effect is generally important 
when the fluctuations are on superhorizon scales, i.e., the characteristic wavelength is larger than 
the Hubble radius. It is small in single-field slow-roll models because the linear approximation is 
valid to high accuracy [20], generically suppressed by the slow-roll parameter r/y defined in Eq. (13). 

For multi-field models, however, the contribution to the energy momentum tensor from some 
of the fields can be highly nonlinear as depicted in Fig. 3. 



B. Non-Gaussianity from multi-field dynamics 
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FIG. 3: An illustration of the energy density configuration in the multi-field case. The density of the 
A- matter/field pa may vary nonlinearly without significantly affecting the total energy density. 

The important property of non-Gaussianity in this case is that it is always of the spatially local 
type. Namely, to second order in nonlinearity, the curvature perturbation will take the form [21], 

ll c (x) = K c , (x) + lf%?Kl (x) , (33) 

where 1Z c ,o is the Gaussian random field and /]y£ al is a constant representing the amplitude of 
non-Gaussianity. The factor 3/5 in front of /Jyff 1 is due to a historical reason. The reason why it 
is of local type is simply causality: No information can propagate over a length scale greater than 
the Hubble horizon scale. 

Observationally, this type of non-Gaussianity can be tested by using the so-called squeezed type 
templates where one of the wavenumbers, say k\ in the bispectrum B(k±, k2, k%) is much smaller 
than the other two, ki <C &2 ~ ks [14], and there are a few observational indications that /yv£ al 
is actually non-vanishing. For example, the WMAP 7 year data analysis gave a one-sigma bound 
11 < f x £f l < 53 (68% CL) [13]. 

IV. 8N FORMALISM 

As mentioned in Sec. II, the 5N formalism is a powerful tool to evaluate the comoving curvature 
perturbation on superhorizon scales. It then turned out that it can be easily extended to the 
evaluation of nonlinear, non-Gaussian curvature perturbations [22, 23]. Let us recapitulate its 
definition and properties: 

(1) SN is the perturbation in the number of e-folds counted backward in time from a fixed final 
time, say t = tf, to some initial time t = t{. 

(2) The final time tf should be chosen such that the evolution of the universe has become unique 
by that time, i.e., the universe has reached the adiabatic limit. Then the hypersurface t = tf 
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should be identified with a comoving (or uniform density) slice, and the initial hypersurface 
t = ti should be identified with a flat slice. 

(3) 6N is equal to the conserved (nonlinear) comoving curvature perturbation on superhorizon 
scales at t > tf. 

(4) By definition, it is nonlocal in time. However, because of its purely geometrical definition, 
it is valid independent of which theory of gravity one considers, provided that the adiabatic 
limit is reached by t = tf. 




cntropy/isocurvature — > adiabatic 



FIG. 4: Three different types of SN. The field space (cf>i, 02) in the figure represents the degrees of freedom 
in the initial condition of the universe. The adiabatic limit is defined to be the stage by which all the 
trajectories converge to a unique one. 

There are various kinds of sources that generate 6N. They may be classified into three types, 
as depicted in Fig. 4. The left one describes a perturbation along the evolutionary trajectory of 
the universe. This case is the same as that of single-field slow-roll inflation, in which the comoving 
curvature perturbation is conserved all the way until it re-enters the horizon. The middle one is 
the case when a small difference in the initial data develops into a substantial difference in 8N. 
Typically this is realized when there is some instability orthogonal to the trajectory, like the case 
when the scalar field moves along a ridge. This type of sources of 5N usually induces a feature 
in the spectrum and/or bispectrum of the curvature perturbation. The right one represents the 
case when the perturbation orthogonal to the trajectory does not contribute to the curvature 
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perturbation until or after the end of inflation, but 5N is generated due to a sudden transition 
that brings the universe into an adiabatic stage. Typical examples are curvaton models [24-26] 
and multi-brid inflation models [27, 28]. 

Here, for the sake of completeness, let us present the precise definition of the nonlinear SN 
formula. See Fig. 5. It is based on the leading order approximation in the spatial gradient expansion 
or the separate universe approach [22], where spatial derivatives are assumed to be negligible in 
comparison with time derivatives. At leading order of the spatial gradient expansion, if we express 



the spatial volume element as \A 3 H = a 3 (t) exp[37£(i, x)] where a(t) is the scale factor of a fiducial 
homogeneous and isotropic universe, we easily find that the perturbation in the number of e-folds 
along a comoving trajectory between two hypersurfaces t = t\ and t = t 2 is given by 



6N(t2,t 1 ;x i ) = n(t 2 ,x i )-n(t 1 ,x i ), 



(34) 



where x l are the comoving coordinates. Here we note that this is purely a geometrical relation. It 
has nothing to do with any equations of motion. 



S G (f ): comoving/uniform density 




p(t 2 ) = const. 

WJ = o 



FIG. 5: Definition of nonlinear SN. It is defined as the perturbation in the number of e-folds from an initial 
flat slice to a final comoving slice. 



First we fix the final hypersurface t = t 2 . It should be taken at the stage when the evolution 
of the universe has become unique. That is, there exists no isocurvature perturbation any longer 
that could develop into an adiabatic perturbation at later epochs. Thus the comoving curvature 
perturbation is conserved at t > t 2 . In the context of the concordance ACDM model of the universe, 
this corresponds to the final radiation-dominated stage of the universe. 

Next we choose the initial slice t = t\. It should be chosen to be flat. Here 'flat' means that 
the perturbation in the spatial volume element vanishes. Namely, the flat slice is defined as a 
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hypersurface on which 1Z = 0. We note that despite its name, the scalar curvature vanishes only 
in the linear theory limit: It is non-vanishing in general in the nonlinear case. 

Applying the above choice of the initial and final hypersurfaces to Eq. (34), it is trivial to see 
that we have 

SN(t 2 ,h-x i ) = n c (t 2 ,x i ). (35) 

Now by assumption 1Z C is conserved at t > t 2 - So it is the quantity we want to evaluate. This 
completes the derivation of the nonlinear 5N formula. 

As mentioned above, since Eq. (34) is a pure geometrical relation, so is the nonlinear SN 
formula (35). This is the reason why it can be applied to any theory of gravity as long as it is a 
geometrical (i.e. general covariant) theory. 

Of course, the above definition tells us nothing about how to evaluate it in practice. In this 
respect, we have a very fortunate situation in the case of inflationary cosmology. It is the fact 
that the evaluation of the quantum fluctuations of the inflaton field, whether it is single- or multi- 
component, can be most easily done in a gauge in which the time slicing is chosen to be flat [12]. 
Thus we can choose the initial slice to be an epoch when the scale of our interest has just exited 
the horizon during inflation. Let the fluctuations of a multi-component scalar field on the flat slice 
at t = t\ to be S(j) a . Then assuming that the values of the scalar field determine the evolution of 
the universe completely, which is the case for slow-roll inflation, the nonlinear 5N can be simply 
evaluated as 

SN = N{4> a + 5<j) a ) - N(0 a ) , (36) 

where N((f) a ) is the e-folding number of the fiducial background. In particular, to second order in 
5(j) a , we obtain 

r)N 1 8 2 N 

n ' = SN =ar sr + 2We¥ SrS * b + -' (37) 

Comparing this with Eq. (33), we see that the curvature perturbation takes a bit more complicated 
form that the simplest form. Nevertheless if we consider the bispectrum, i.e., the Fourier compo- 
nent of the three-point function (1Z c (xi)lZ c (x2)1Z c (x3)) , we find there is a quantity that exactly 
corresponds to ffif 1 defined in Eq. (33). Namely [23], 

3 f iocai _ G ab G cd N a N bc N d N= dN_ d 2 N 

5 Jnl ~ 2(||ViV||2) 2 ' a ~ 84> a b ~ drdcfP' 1 ' 

Before concluding this section, we mention the fact that the SN formalism does not require the 
scalar field fluctuations to be Gaussian. In fact, except for the last equation in the above, Eq. (38) 
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which assumes the Gaussianity of 5(j) a , the general SN formula (36) or its second order version (37) 
can be used for non-Gaussian 5cf) a [29]. Such a case may happen, for example, in multi-field DBI 
inflation. 

V. SUMMARY 

It has been about 30 years since the inflationary universe was first proposed, and there is increas- 
ing observational evidence that inflation did take place in the very early universe. Among others, 
the measured CMB temperature anisotropy is fully consistent with the predictions of inflation 
that the primordial curvature perturbation spectrum is almost scale-invariant and it is statistically 
Gaussian. 

Inflation also predicts a scale-invariant tensor spectrum, and if the energy scale of inflation is 
high enough as in the case of chaotic inflation, the tensor-scalar ratio r can be as large as 0.1. If 
this is the case, the tensor perturbation will be detected in the near future, and it will confirm not 
only the inflationary universe but also quantum gravity. 

Even if the tensor perturbation will not be detected, there may be other interesting signatures 
of inflation. Non-Gaussianity from inflation is attracting attention as one of those signatures that 
can distinguish or constrain models of inflation significantly. 

We discussed that the origins of primordial non-Gaussianities may be classified into three cat- 
egories, according to different length scales on which different mechanisms are effective: 

(1) Quantum theoretical origin on subhorizon scales during inflation. 

(2) Classical nonlinear scalar field dynamics on superhorizon scales during or after inflation. 

(3) Nonlinear gravitational dynamics after the horizon re-entry. 

In particular we argued that non-Gaussianities in the second case are always of spatially local 
type. We then mentioned that there are three different kinds of situations in which such local non- 
Gaussianities can be generated, and described in some detail a very efficient method to compute 
them, namely, the 5N formalism. 

Apparently identifying properties of primordial non-Gaussianities in the observational data is 
extremely important for understanding the physics of the early universe. Here we mentioned only 
the bispectrum or the 3-point function. But if it is detected, higher order n-point functions may 
become important as a model discriminator. Other types of non-Gaussianity discriminators may 
also become necessary. 
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What is important is that we are now beginning to test observationally the physics of the very 
early universe, the physics at an energy scale closer to the Planck scale, at a scale that can never 
be attained in high energy accelerator experiments. 

Cosmology has become not only a precision science, but now it constitutes a truly indispensable 
part of fundamental physics. General relativity is the backbone of cosmology. I wonder what 
Einstein would say if he were here in this very exciting era - 100 years after he visited Prague. 
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